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The spatial search problem on regular lattice structures in integer number of dimensions d > 2 
has been studied extensively, using both coined and coinless quantum walks. The relativistic Dirac 
operator has been a crucial ingredient in these studies. Here we investigate the spatial search 
problem on fractals of non-integer dimensions. Although the Dirac operator cannot be defined on 
a fractal, we construct the quantum walk on a fractal using the flip-flop operator that incorporates 
a Klein-Gordon mode. We find that the scaling behavior of the spatial search is determined by 
the spectral (and not the fractal) dimension. Our numerical results have been obtained on the 
well-known Sierpinski gaskets in two and three dimensions. 

PACS numbers: 03.67.Ac 
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INTRODUCTION 



The spatial search problem is to find a marked object 
from an unsorted database of size N spread over distinct 
locations, with the restriction that one can proceed from 
any location to only its neighbors while inspecting the 
objects. Classical algorithms for the unsorted database 
search can do no better than inspect one location after 
another, and so are 0{N). On the other hand, quantum 
algorithms can do better by working with a superposi- 
tion of states. The quantum spatial search problem has 
been studied extensively in recent years, in a variety of 
geometries ranging from a single hypercube to regular 
lattices in various integer dimensions (see for example, 
Refs. I1M3). These investigations used local translation- 
ally invariant quantum walk operators [sl, Q , arising from 
relativistic quantum mechanics in some form, to obtain 
their results. Here we investigate the spatial search on 
a fractal lattice, which has neither a translational sym- 
metry nor a straightforward structure to implement the 
relativistic Dirac operator. The non- integral value of the 
fractal dimension also lets us explore the dependence of 
the spatial search on the geometry of space in more gen- 
erality. 

Our algorithmic strategy for the quantum search is to 
construct a Hamiltonian evolution, whereby the kinetic 
term diffuses the amplitude distribution all over the space 
and the potential term attracts the amplitude distribu- 
tion toward the marked vertex [10] . In its discrete form, 
the kinetic term is realized as the quantum walk oper- 
ator W, explicitly defined in Section III. A relativistic 
walk provides the fastest diffusion, and is expected to 
produce the quickest search. The best potential term 
is the one that provides maximum contrast between the 
marked vertex and the rest. In its discrete form, it is the 
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where r is the time-step size. The algorithm alternates 
between the oracle and the walk operators, yielding the 
time evolution 



mx;h,t2)) = [W'^Rf^i^{x;0,0)) . 



(2) 



Here ^2 is the number of oracle calls, and ti is the num- 
ber of walk steps between the oracle calls. Both have 
to be optimized, depending on the spatial distribution 
of the vertices, to concentrate the amplitude distribution 
toward the marked vertex as quickly as possible and to 
solve the spatial search problem. 

The iterative evolution of Eq.(|2]) redistributes the am- 
plitude at each vertex over itself and its neighbors at 
every step. Unitarity of the evolution means that the 
eigenvalues have unit magnitude, and so the results of 
the algorithm are periodic in time. Grover's algorithm 
is the special case where one can move from any vertex 
to any other vertex in just one step, which can be inter- 
preted as the d ^ oo limit of spatial search [6] . Unlike 
Grover's algorithm, the maximum probability of being at 
the marked vertex, P, does not reach the value 1 for a 
generic spatial search. Augmenting the algorithm by the 
amplitude amplification procedure [11|, the marked ver- 
tex can be found with probability 6(1), and the overall 
complexity of the algorithm is then characterized by the 
effective number of oracle calls t2/VP. 

We have argued [6] that the spatial search in d dimen- 
sions obeys two lower bounds. 



t2 > max{(i7Vi/^,7rVA^/4} . 



(3) 
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They, respectively, follow from the distinct physical prin- 
ciples of special relativity (i.e. the walk must be able to 
travel across the lattice) and unitarity (i.e. Grover's algo- 
rithm is optimal in the absence of any constraint on walk 



movements). The best spatial search algorithms, there- 
fore, appear in the framework of relativistic quantum me- 
chanics. Moreover, an analogy with statistical mechanics 
of critical phenomena in different dimensions, having an 
interplay of multiple dynamical features, provides an un- 
derstanding of the scaling behavior of the spatial search. 
For (i > 2, the latter bound wins, and the spatial search 
is 0{\fN) differing only in scaling prefactors. The two 
bounds cross in the critical dimension d = 2, where loga- 
rithmic corrections to scaling behavior are expected and 
have indeed been observed [2|-[5|,[3- For d < 2, the former 
bound dominates and should govern the scaling behavior 
of the spatial search. But the relativistic evolution op- 
erator is power-law infrared divergent there (as J dfik/k'^ 
in the continuum formulation), which may modify the 
scahng behavior compared to the bound. 

The preceding arguments do not provide any new in- 
sight for (i = 1, where both classical and quantum spatial 
search are 0{N). We need non- integral values of d to 
convincingly test the analogy of the spatial search with 
critical phenomena in statistical mechanics, and to that 
end we have investigated the spatial search on fractal lat- 
tices. Many candidates for fractal lattices are available 
with 1 < d < 2, and for them we expect the spatial search 
to have a power-law scaling form, say 0{N^). We want 
to determine how the exponent s depends on d, in par- 
ticular whether it indeed equals l/(i, and we specifically 
work with fractals known as Sierpinski gaskets. 



II. QUANTUM WALK 
IN A FRACTAL GEOMETRY 



Fractals are self-similar structures, and are conve- 
niently defined using a recursive scheme in an embedding 
space, say of Euclidean dimensions dE > d. (This is in 
contrast to regular lattices generated by repetitive unit 
cells in translationally invariant systems, where ds = d.) 
For Sierpinski gaskets, the recursive structure is a regular 
simplex. That is an equilateral triangle for d^; = 2, and 
Fig. 1 illustrates the recursive generation of the Sierpinski 
gasket as a function of the stage number S. 

At stage S, the linear extent of Sierpinski gaskets is 
L = 2^ . Since a simplex has d^; + 1 vertices, the total 
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number of vertices N in Sierpinski gaskets is given by, 

(d£; + l)((d£; + 1)^ + 1) 
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(4) 



Thus the fractal (also called the Hausdorff ) dimension of 
Sierpinski gaskets is 



d = \n{dE^l)/\n2 . 



(5) 



Dynamical phenomena on fractals have been investi- 
gated before [l2|, [HI , in particular the classical random 
walk and the vibrational problem. Instead of the trans- 
lational symmetry of regular lattices, self-similar fractals 
have the dilation (or the scale) symmetry. That leads 
to power-law scaling behavior for physical observables, 
and the critical exponents for various well-known frac- 
tals have been calculated. The most common operator 
appearing in the dynamical equations is the Laplacian, 
and the critical behavior is then governed by the scaling 
of its spectral density. 
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This result can be understood as the reciprocal space for 
the Laplacian (or the Brillouin zone) having an effective 
dimension ^5, which is called the spectral (or the fracton) 
dimension |12|, [l3|. In the case of Sierpinski gaskets. 
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FIG. 1: Sierpinski gasket generation at various stages. 



which cannot exceed 2 although d can. 

For regular lattices, dE = d = ds, whereas it is ob- 
served that dE ^ d > ds in general. Note that both the 
position and the reciprocal space describe the same num- 
ber of degrees of freedom in the same number of embed- 
ding dimensions. So d ^ ds means that the arrangement 
of vertices is different (in terms of connectivity as well 
as spacings) in the position and the reciprocal spaces. 
For the Sierpinski gaskets in two and three embedding 
dimensions, the relevant dimensions are: 

dE = 2 : d= 1.5849625 . . . , 4 = 1.3652124 . . . , 
dE = 3 : d = 2, 4 = 1-5474112.... (8) 

Both the classical random walk and the vibrational 
problem involve the same Laplacian in their spatial de- 
pendence, although they differ in their temporal depen- 
dence (first and second derivative respectively). So their 
critical exponents depend on the same ds, but the de- 
pendence differs in phenomena involving temporal cor- 
relations. For the relativistic quantum walk problem we 
are addressing, the appropriate scaling behavior corre- 
sponds to that for the vibrational problem (and not the 
classical random walk). 

Now quantum walk diffuses at relativistic speeds on a 
hypercubic lattice when the Dirac operator is used for 
time evolution. But it is unclear how to extend the rele- 
vant Clifford algebra to non-integer dimensions [MJ. So 



our earlier methodology for the spatial search on hyper- 
cubic lattices [sl^SH] cannot be implemented on fractals, 
and we need a different implementation of a quantum rel- 
ativistic walk. It is known that the flip-flop walk diffuses 
relativistically [2], and can be easily applied to graphs 
where all the vertices have the same connectivity. That 
happens to be the case for Sierpinski gaskets, which have 
2dE neighbors for every vertex. We therefore first con- 
struct the flip-flop walk on Sierpinski gaskets, in terms 
of independent walk steps along every link, and then use 
it to study the spatial search problem. 



III. THE FLIP-FLOP WALK 

The flip-flop walk [2| is constructed in the joint Hilbert 
space containing both the position and the link degrees 
of freedom, 



IV^(f,0) = ^a^,z|f)0|O enN^nk=% 



search 



(9) 



Every link attached to a given vertex is assigned its own 
direction pointing away from the vertex, so the dimen- 
sion k of Hk is the number of neighbors for the vertex. 
(For ease of implementation, we take this number to be 
the same for each vertex.) The walk step consists of a 
shift operation followed by a mixing of the link ampli- 
tudes, i.e. W = GS. The shift operation first propagates 
the quantum amplitude along its link direction, and then 
reverses the link direction. 



\x)^\l)^\x^l)^\-l) 



(10) 



Note that when the graph has a link / at the vertex x, 
it automatically has a link — / at the vertex x -\- 1. The 
mixing of links at every vertex is the reflection in the 
mean operation (also referred to as the Grover coin [9]), 
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Both S and G are reflection operators, S^ = 1 = G^^ 
and the walk spreads on the lattice because these two 
operators do not commute. S is block diagonal on the 
links, and its eigenstates are the combinations |x, /) =b |x + 
/, — /) for every link, with eigenvalues ±1. Explicitly, 

.5=0 [{\xJ) + \x+i,-iMxJ\ + {x+i,-i\)-l] . (12) 

links 

On the other hand, G is block diagonal on the vertices. 



G= \-AT.\-^h){T.(-^i'\)-i 



The unbiased uniform superposition state, 

x,l 



(14) 



is invariant under both S and G (i.e. both operators 
have eigenvalue +1). Furthermore, the oracle R trivially 
commutes with G. 

When the lattice has translational invar iance, S can 
be simplified by Fourier transforming to the momentum 
space, and then the eigenspectrum of the walk can be 
obtained by diagonalizing a k x k matrix describing the 
operator W for each value of the momentum [2]. Such 
an analysis, however, is not possible for a fractal. 



A. Flip-flop walk on Sierpinski gaskets 

We first implement the flip-flop walk on the Sierpinski 
gasket in two dimensions. As can be seen from Fig. 1, all 
the internal vertices of the gasket have four neighbors, 
whereas the three corner vertices have only two neigh- 
bors. We increase the number of neighbors for the corner 
vertices to four, by providing a periodic wraparound as 
illustrated in Fig. 2. Furthermore, to conveniently label 
the vertices and to keep track of their neighbors, we em- 
bed the gasket in a two-dimensional rectangular grid as 
shown in Fig. 3. 

In this configuration, there are six possible walk direc- 
tions in general, as depicted in Fig. 4, out of which only 
four are valid for any given vertex. We separate the ver- 
tices in to six types, also shown in Fig. 4, depending on 
which of the four walk directions are applicable to them. 
Note that VO, VI and V2 are internal vertices, whereas 
V3, V4 and V5 are the corner vertices. 

With this infrastructure, we can now explicitly write 
down the operations involved in the flip-flop walk. The 



(13) 




For every vertex, it has one eigenvalue +1 corresponding 
to the isotropic distribution, and k — 1 eigenvalues —1 
corresponding to other orthogonal distributions. 



FIG. 2: Wraparound for the corner vertices of the two- 
dimensional Sierpinski gasket. 
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FIG. 3: Embedding of the two-dimensional Sierpinski gasket 
in a rectangular grid. 



shift operation is: 
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The mixing of link amplitudes occurs only among the 
four valid directions at each vertex. It can be described 
by the operator matrix, 
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These definitions can be easily extended for fiip-fiop 
walks on Sierpinski gaskets in any number of dimensions, 
with k = 2dE' 



B. Flip-flop walk as a relativist ic propagator 

To understand the time evolution provided by the fiip- 
fiop walk operator, we describe it for a hyper cubic lattice 
in an integer number of dimensions, in a notation that 
identifies Hk with internal degrees of freedom. 

For (i = 1, the shift operator is 



S = ^(l^ai) 



^ hx+l)(x| (g)Cr_ + |x-l)(x|(g)Cr+ , (17) 



and the link- mixing operator is G^^^^^ = 1 (g) <Ji. To- 
gether they produce a directed walk, where the two spinor 
components propagate independently and at full speed in 
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FIG. 4: Above: possible walk directions, and below: various 
vertex types, for the two-dimensional Sierpinski gasket. 



opposite directions. 
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(18) 
This is a representation of the massless Dirac propagator, 
in the basis where the two spinor components describe 
the two chiral degrees of freedom in d = 1. It is also the 
representation of the massless scalar Klein-Gordon prop- 
agator, where the two parts of the d'Alembert solution 
travel independently in opposite directions. 

The situation in higher dimensions is more non-trivial. 
For a hypercubic lattice in d dimensions, size L = 7V^/^ 
and connectivity k = 2d. The matrix W = GS is uni- 
tary, with eigenvalues of unit magnitude and orthogonal 
eigenvectors. Its spectrum can be obtained as follows. 
The momentum components on a periodic lattice are 
ki = 27r/^/L, with li G {0,1,..., L — 1}. Select a spe- 
cific set of momentum components {ki}. Then pick any 
two-dimensional subspace, from the {d — 1) linearly in- 
dependent choices (xi,X2), (xi,X3), . . . , (xi,^^^). In this 
subspace. 



(19) 



and Gij = —5ij + 1/d. Corresponding to each di- 
rection, S has eigenvalues A = ±1 with eigenvectors 
(l,Ae*^^). From the degenerate pair of eigenvectors in 
every two-dimensional subspace, one can construct com- 
binations that are orthogonal to (1,1,1,1). Acting on 
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these combinations, G reduces to — /, and so the com- 
binations are eigenvectors of W with eigenvalues —A. 
Exphcitly, they are (1 + \e'^\\e'^-{l + Ae^^^),-(1 + 
Ae*^''), — Ae*^^(l + Ae*^^)), with vanishing components 
outside the two-dimensional subspace. The {d — 1) lin- 
early independent choices of two-dimensional subspaces 
thus yield {d — l)-fold degenerate eigenvalues ±1 for W . 
The remaining two eigenvalues of W are a complex conju- 
gate pair e'^'^^^^^\ and Tr(W) determines their real part 
as co^{uj(ki)) = \Y.i cos(/ci). 

We observe that the ±1 eigenvalues of the flip-flop walk 
correspond to local modes that can be restricted to two- 
dimensional subspaces, whereas the complex conjugate 
pair of eigenvalues provides relativistic propagation with 
speed 1/yd. With this break up of the evolution modes, 
we can surmise that the flip-flop walk contains a sin- 
gle first-order form [15] of the complex massless Klein- 
Gordon operator, which makes its dynamics relativis- 
tic. The scalar Klein-Gordon operator can be mathe- 
matically constructed in any number of dimensions, even 
non- integer ones, unlike the spinor Dirac operator. Our 
simulation results in the next section demonstrate that 
the flip-flop walk indeed achieves the desired goals for 
the spatial search problem, and it would be interesting 
to explore its usefulness in other relativistic problems in- 
volving scalar fields in non-integer dimensions (or curved 
spaces). 



IV. SPATIAL SEARCH 
USING THE FLIP-FLOP WALK 

We first simulated the spatial search algorithm, for a 
single marked vertex on Sierpinski gaskets, as defined by 
Eq.([2]). We chose the initial state to be the unbiased uni- 
form superposition state, as is customary. Our numerical 
results are described later in the section. In particular, 
we needed ti = 2 for the search to succeed with reason- 
able probability (ti = 1 did not work as well, and ti = 3 
was far worse), and both the success probability and the 
number of oracle calls showed power-law dependence on 
the database size N. Note that our spatial search al- 
gorithm is different from that of Ref. [2] for two reasons 
(although we use the same flip-flop walk): (i) we allow for 
values of ti 7^ 1, and (ii) our oracle flips only the sign of 
the marked vertex amplitude, whereas Ref. [2] applies the 
link-mixing operator G^^^ at the marked vertex together 
with the sign flip. 

Since the flip-flop walk corresponds to a massless evo- 
lution propagator, this algorithm suffers from infrared 
divergence and is not optimal. Tulsi showed how to 
eliminate the infrared divergence of the spatial search 
algorithm by controlling the evolution operators using 
an ancilla qubit [4]. The ancilla control traps the quan- 
tum walk at the marked vertex and enhances the algo- 
rithm's success probability. Explicitly, in the joint space 
'Hanciiia ® '^search, the |1) (g) |x, /) statcs are the original 
search space states, the |0) ^ \x = t/l) state is the new 
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FIG. 5: Logic circuit diagram for Tulsi's controlled quantum 
spatial search algorithm. R and W are the binary oracle and 
the quantum walk operator respectively. We use the general- 
ization with two walk steps after each oracle. 



trap state that develops a nonzero amplitude, and the 
|0) |x 7^ t, /) states maintain zero amplitudes. This fea- 
ture can be interpreted as the introduction of an effective 
mass at the marked vertex in the Hamiltonian Q, or as 
the introduction of a self- loop at the marked vertex in the 
graph [16]. Tulsi's scheme is illustrated in Fig. 5, where 
the ancilla operators are 
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(20) 



and the algorithm evolves the initial state |1) (g) \s) to the 
target state |^) \tj) with \S) = X]|l). Note that for 
the generic state corresponding to Eq.®, the probability 



of being at the marked vertex is Xl^d^^- 



(0) 



--t,i\ 



.(1) 



=tJ 



n 



where the superscript refers to the ancilla value. 

For J = 0, Tulsi's algorithm reduces to spatial search 
described by Eq.([2]). It finds the marked vertex with 
probability Pq = e(7V"^) using Qo = 0{N^) oracle calls. 
Tulsi showed that, with cos 5 = Q{^/Po)^ the algorithm 
increases the probability of finding the marked vertex to 
Ps = 0(1) without changing the scaling of oracle calls 
Qs = Q{N^) [4]. More explicitly, the algorithm largely 
confines the evolution of the quantum state to a two- 
dimensional subspace of Hat, whereby 
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Here B = Bq = Q{N^^'^) is a second moment constructed 
from the eigenspectrum of W, and characterizes the in- 
frared divergence of the problem. Also, Eq. (|2T]) relates 
the scahng exponents. 



a = 2b-l . 



(23) 



The optimal value of the ancilla control parameter is 
obtained by minimizing the algorithmic complexity. 
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FIG. 6: (Color online) A snapshot of the probability distribu- 
tion for the spatial search on the two-dimensional Sierpinski 
gasket, when the marked vertex attains its peak probability. 
The data are for S — 6 with ancilla control. 
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FIG. 7: (Color online) Time evolution of the marked vertex 
probability for the spatial search on the two-dimensional Sier- 
pinski gasket. The data are for S = 6. Those without ancilla 
control are shown as the dotted curve, and those with ancilla 
control are shown as the continuous curve. 
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In what follows, we check these expectations regarding 
the scaling behavior of the spatial search, for Sierpinski 
gaskets in two and three dimensions. 



A. Simulation results for ds = '2 

We simulated the spatial search algorithm on the two- 
dimensional Sierpinski gasket, for stage 5* = 4 — 14 cor- 
responding to A/" = 123 — 7174455. To keep the memory 
requirements under control, we kept track of the irregular 
pattern of vertices and their neighbors using hash tables. 
We first checked that the algorithm did produce an ap- 
proximately periodic evolution for the probability at the 
marked vertex, and the probability peak was reasonably 
sharp for the parameters we chose. For S = 6 {N = 1095) 
and the marked vertex in the center of the gasket, a snap 
shot of the probability distribution is presented in Fig. 6, 
and the probability evolution at the marked vertex is 
shown in Fig. 7. We find that the probability evolution 
is smoother with ancilla control than without it, imply- 
ing that ancilla control improves the confinement of the 
evolution of the quantum state to a two-dimensional sub- 
space oH-Ln- 

Implementing the spatial search algorithm without an- 
cilla control was quite straightforward. But to implement 
Tulsi's algorithm, we needed the optimal value of cos (5, 
and we did not have any analytical estimate for it. So we 
just assumed Eqs. (f2T|) and ([24|) to be valid, and extracted 
(cos5)opt from the result for Pq at each stage S. That 
worked well, and we obtained P^ ^ 0.5 as predicted. 

Lack of translational invariance for a fractal means 
that the spatial search results depend on the location of 
the marked vertex. In the absence of ancilla control, each 
of Qo, ^0 and Qo/V^ showed variation up to a factor 
of 2.2, as we moved around the position of the marked 



vertex in the gasket. With ancilla control, variation in 
Q5 was somewhat smaller, up to a factor of 1.6, whereas 
Ps ~ 0.5 remained valid (within 10%) everywhere as ex- 
pected. We also observed that the largest variation in 
the spatial search results was between the marked ver- 
tex positions at the corner and the center of the gasket. 
More importantly, when we extracted the scaling behav- 
ior from our results, holding the relative position of the 
marked vertex in the gasket (i.e. x/L) fixed, we found 
that the variation was essentially in the prefactors and 
the scaling exponents changed by less than a percent. 
Henceforth we opted to keep the marked vertex at the 
center of the gasket, and all our results that follow cor- 
respond to that choice. 

Our data for the number of oracle calls Q and the 
corresponding peak probability P are presented in Figs. 8 
and 9, both without and with ancilla control. The figures 
also show simple scaling fits that were performed in the 
range 5 = 8 — 14. 

In the absence of ancilla control, our fits give: 

logsQo = -1.066 + 0.730 log2A^ (err = 0.004), 
log2 Po = 1.995 - 0.440 logs N [err = 0.030).(26) 

Here ^^err^^ refers to the r.m.s. deviation in the data from 
the fit. The dominant error in the asymptotic scaling 
parameters is the systematic error due to finite values of 
N. So we quote as our error estimates, the difference in 
the numbers for fits corresponding to 5 = 7 — 14 and 
5 = 8 - 14. Then the results of Eq.([26l) translate to the 
scaling behavior: 



Qo = 0.478(13) Ar0-730(2) ^ 
Po = 3.99(20) Ar-o-440(4) 



(27) 



The scaling exponents satisfy the relation Eq. ([23|) well, 
which supports Tulsi's analysis criterion that the quan- 
tum state evolution is largely confined to a two- 
dimensional subspace of Hn- 




log2(N) 

FIG. 8: (Color online) Scaling of the number of oracle calls for 
the spatial search on the two-dimensional Sierpinski gasket. 
The open and filled symbols denote data without and with 
ancilla control respectively. The linear fits are for the data 
from S = 8to S=U. 



constant close to 0.5. Both these features fully agree with 
Tulsi's predictions. We estimate the scaling behavior as: 

Qs = 0.599(8) 7V0-724(i) ^ p^ ^ 0.4903(2) . (29) 

Also, direct fits for the complexity behavior give: 

log2(Q5/\/^) = -0.334 + 0.729 logsA^ (err = 0.015), 
Q6/^/P6 = 0.794(13) 7V0-729(i) . (39) 

We observe that the exponent b ^ j-^ = 0.73248676 . . .. 
That clearly obeys the bound in Eq.(|3]), since l/d < l/dg. 
But more than that, our results imply that, for a fractal, 
the relevant length scale for the spatial search is not its 
linear extent L ^ A^^/^, rather it is the inverse of the 
linear extent of its reciprocal space N'^/^^ . Thus it is the 
smaller spectral dimension which governs the dynamics 
of the search process on fractals. 
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FIG. 9: (Color online) Scaling of the peak probability for 
the spatial search on the two-dimensional Sierpinski gasket, 
without (top) and with (bottom) ancilla control. The fits 
are for the data from 5* = 8 to 5* = 14, linear (top) and 
exponential (bottom) respectively. 



The fits for our results with ancilla control, with cos S 
determined from the value of Pq, give: 

logsQ^ = -0.739 + 0.724 log2 TV (err = 0.017), (28) 
P5 = 0.4903 + 0.0471 e"°-°°^^^^ (err = 0.0023). 

As displayed in Fig. 8, Qs is only marginally higher than 
Qo, and the scaling exponent b is essentially unchanged. 
On the other hand. Fig. 9 shows that Ps approaches a 



B. Simulation results for ds = 3 

To reinforce our observations, we extended our sim- 
ulations of the spatial search algorithm to the three- 
dimensional Sierpinski gasket. In this case, the recursive 
structure is a tetrahedron, whereas the linear extent at 
stage S remains L = 2^ . Each internal vertex of the 
gasket now has six neighbors, and every corner vertex 
has a periodic wraparound with the other three corner 
vertices. Furthermore, there are twelve possible walk di- 
rections in general, out of which only six are valid at any 
given vertex. In our explicit operations, we used six types 
of internal vertices, four types of corner vertices, and the 
link- mixing operator G(^=^) . 

We obtained results for stage 5* = 3 — 10 corresponding 
to N = 130 — 2097154, and the marked vertex at the 
center of the gasket. Although not as extensive as the 
results for the (Ie = 2 case, these are sufficient to discern 
the patterns in the scaling laws. 

We again first simulated the algorithm without ancilla 
control, and then used (cos (5)opt extracted from the value 
of Po for each stage S to simulate the algorithm with 
ancilla control. Our results are plotted in Figs. 10 and 11, 
together with the simple scaling fits that were performed 
in the range 5* = 5 — 10. Without ancilla control, our fits 
give: 

logsQo = -0.767 + 0.641 logsA^ (err = 0.002), 
log2 Po = 1.39 - 0.265 log2 N {err = 0.036). (31) 

Again, ^^err'''' refers to the r.m.s. deviation in the data 
from the fit. For the error estimates on the fit parame- 
ters, we quote the difference in the numbers for fits cor- 
responding to 5 = 5 — 10 and 5 = 6 — 10. Then the 
corresponding estimates of the scaling behavior are: 



Qo = 0.588(1) Ar0-64i(i) ^ 
Po = 2.62(25) Ar-0-265(T) 



(32) 




0.518 {err = 0.001). 



(33) 



log2(N) 

FIG. 10: (Color online) Scaling of the number of oracle calls 
for the spatial search on the three-dimensional Sierpinski gas- 
ket. The open and filled symbols denote data without and 
with ancilla control respectively. The linear fits are for the 
data from S = 5 to S = 10. 
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FIG. 11: (Color online) Scaling of the peak probability for 
the spatial search on the three-dimensional Sierpinski gasket, 
without (top) and with (bottom) ancilla control. The fits are 
for the data from 6' = 5tOAS' = 10, linear (top) and constant 
(bottom) respectively. 



Once more, the scaling exponents reasonably satisfy the 



relation Eq.(|23l), and 6 ?^ ^ = 0.64624063 . . .. 
With ancilla control, our fits yield: 

log2Q<5 = -0.817 + 0.647 logsA^ (err = 0.003), 



We observe that Qs is only slightly larger than Qo, the 
scaling exponent b is essentially the same, and Ps is a 
constant close to 0.5, all consistent with Tulsi's predic- 
tions. We estimate the scaling behavior as: 

Qs = 0.568(4) Ar0-647(i) ^ p^ ^ 0.518(1) . (34) 

Moreover, direct fits for the complexity behavior give: 



-0.343 + 0.647 log2 N {err = 0.002), 

0.788(3) 7V0-647(i) (35) 






Here the fact that 6 ~ ^ is an important finding. 
In this case, the embedding dimension d^; = 3 and the 
Hausdorff dimension d = 2. Were either of them control- 
ling the scaling behavior of the spatial search, we would 
have found Q = 0(>/7V). Our clearly distinct result once 
again tells us that the dynamics of the spatial search is 
governed by the spectral dimension. 



V. CONCLUSION 

We have described how the flip-flop walk contains 
within it a relativistic Klein-Gordon propagator, and 
then used it to solve the spatial search problem on frac- 
tals. The analogy with critical phenomena in statistical 
mechanics suggests that, for 1 < d < 2, physical observ- 
ables have power-law scaling behavior. Also, for d < 2, 
the fully relativistic spatial search algorithm suffers from 
infrared divergence, which we have suppressed by Tulsi's 
ancilla controlled version that is equivalent to the intro- 
duction of an effective mass or a self-loop at the marked 
vertex. 

Specifically, we have carried out numerical simulations 
for Sierpinski gaskets in two and three dimensions, keep- 
ing the relative position of the marked vertex in the gas- 
ket fixed, and our results support all the expectations. 
An important finding is that we are able to reach the 
lower bound of A^^/^ oracle calls, where d is the spec- 
tral dimension and not the fractal dimension. These two 
dimensions are not equal for fractals, and we find that 
the smaller spectral dimension governs the scaling ex- 
ponents. Note that the solution of the quantum spatial 
search problem has to be simultaneous in both the po- 
sition space and the reciprocal momentum space. The 
position space evolution proceeds from a uniform distri- 
bution toward a ^-function, and the momentum space 
evolution proceeds from a (5-function toward a uniform 
distribution. Thus it is fully understandable that the 
stricter of the lower bounds (i.e. the smaller value of dg) 
governs the scaling behavior of the problem. 



[1] N. Shenvi, J. Kempe a nd K.B. Whaley Phys. Rev. A67 

(2003) 052307, arXiv:qTl i[^-ph /0210064| 
[2] A. Ambainis, J. Kempe and A. Rivosh, Proceedings 

of the Sixteenth Annual ACM-SIAM Symposium on 

Discrete Algorithms (SODA'05), (SIAM, Philadelphia, 

2005), p. 1099, arXiv:quant- ph/0402107 
[3] A. Patel, K.S. Raghunathan and P. Rungta, Pro- 
ceedings of the Workshop on Quantum Information, 

Computation and Communication (QICC-2005), IIT 

Kharagpur, India (Allied Publishers, 2006), p. 41, 

'arXiv:quant-ph/0506221 
[4] A. Tulsi, Phys. Rev. A78 (2008) 012310, arXiv:0801.0497; 

see also Ph.D. thesis, Indian Institute of Science (2009). 
[5] G. Abal, R. Donangelo, F.L. Marquezino and R. 

Portugal, Math . Struct. Comp. Sci. 20 (2010) 999, 

tarXiv:1001 . 113 91 
[6] A. Patel and M.A. Rahaman, Phys. Rev. A82 (2010) 

032330, arXiv: 1003.0065, 
[7] A. Patel, K.S. Raghunathan and M.A. Rahaman, Phys. 

Rev. A82 (2010) 032331, ar Xiv: 1003 .5564 
[8] Y. Aharonov, L. Davidovich and N. Zagury, Phys. Rev. 

A48 (1993) 1687. 
[9] J. Kempe, Contemporary Physics 44 (2003) 307, 

arXiv:quant-ph/0303081 
[10] L.K. Grover, Pra mana 56 (2001) 333, 

ar Xiv:quant-ph/0109116| 



[11] G. Brassard, P. H0yer, M. Mosca and A. Tapp, in Quan- 
tum Computation and Information, AMS Contempo- 
rary Mathematics Series Vol. 305, eds. S.J. Lomonaco 
and H.E. Brandt (AMS, Providence, 2002), p. 53, 
arXiv:quant-ph/0005055 , 

[12] S. Alexander and R. Orbach, J. Phys. Lett. 43 (1982) 
625. 

[13] R. Rammal and G. Toulouse, J. Phys. Lett. 44 (1983) 13. 

[14] Note that the Clifford algebra cannot be formulated in 
curved spaces either. Inclusion of fermions in general rela- 
tivity needs vielbeins, with one index in the curved space 
and the other in the flat tangent space. The Dirac op- 
erator structure is defined in the fiat tangent space, and 
is contracted with curved space variables using the viel- 
beins. 

[15] A second-order differential equation can always be rewrit- 
ten as two coupled first-order differential equations, e.g. 



ip -\- ip = and 



d Nj 



are equivalent 



_ 1> 

dt v^y ~ V-i o; \^i;j 

descriptions of a scalar wave equation. 
[16] H. Krovi, F. Magniez, M. Ozols and J. Roland, Pro- 
ceedings of the 37th International Colloquium on Au- 
tomata, Languages and Programming (ICALP 2010), 
Lecture Notes in Computer Science Vol. 6198 (Springer, 
Berlin, 2010), p. 540, arXiv:1002.2419. 



